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Abstract 



We consider periodic matrix-valued Jacobi operators. The spectrum of this operator 
is absolutely continuous and consists of intervals separated by gaps. We define the 
Lyapunov function, which is analytic on an associated Riemann surface. On each sheet 
the Lyapunov function has the standard properties of the Lyapunov function for the 
scalar case. We show that this function has (real or complex) branch points, which we 
call resonances. We prove that there exist two types of gaps: i) stable gaps, i.e., the 
endpoints are periodic and anti-periodic eigenvalues, ii) unstable (resonance) gaps, i.e., 
the endpoints are resonances (real branch points). We show that some spectral data 
determine the spectrum (counting multiplicity) of the Jacobi operator. 

1 Introduction and main results 

Consider self-adjoint matrix-valued Jacobi operators J acting on £ 2 (Z) m and are given by 
(Jy)n = anVn+i+Kyn + an-iVn-i, n G Z, y n <E C m , y = (y n ) n& E f (Z) m , m ^ I, (1.1) 

where a n , b n = are p-periodic sequences of the mxm real matrices and det a n ^ for all 
n E Z. It is well known that the spectrum a(J) of J is absolutely continuous and consists 
of non-degenerated intervals [A^_ 1; A~], X^-i < ^ A+,n = 1,...,N < oo. These intervals 
are separated by the gaps 7„ = (A~,A+) with the length > 0. Introduce the fundamental 
mxm matrix-valued solutions (p = (^(^jnez^ = {&n.( z ))nez °f the equation 

a n yn+i + KVn + oj-il/n-i = z Vn, (f Q = -&i = 0, ip x = O = I m , {z, n) G C x Z, (1.2) 

where I m is the identity mxm matrix. We define the monodromy 2m x 2m matrix ^ p and 
the trace T n , n ^ 1 by 
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Introduce the modified monodromy matrix M and the determinant D by 

M = P ^ P P \ P = ( q° j m J , D(z, r) = det(M(z) - rl 2m ), r,zeC. (1.4) 

Let Tj = Tj(z),j G N 2m be the eigenvalues of M{z), where N m = {1, , ..,m}. An eigenvalue 
of M{z) is called a multiplier. 

Theorem (Lyapunov-Poincare) i) The following identities hold 

M T JM = J, J=( y _° Im (1-5) 

D(z,t) = T 2m D{z,T~ 1 ), all z,reC,T^0. (1.6) 

ii) If t(z) is a multiplier of multiplicity d ^ 1 for some z G C (or z G M.), then t~ 1 (z) 
(or t(z)) is a multiplier of multiplicity d. Moreover, each M(z),z G C. has exactly 2m 
multipliers rf 1 (z),j G N m and a (J) = Uf =1 {z G C : \rj(z)\ = 1}. 
Hi) If t(z) is a simple multiplier for some z G C and \r(z)\ = 1, then r'(z) ^ 0. 

It is well known that the spectrum of the scalar (m = 1) Jacobi operator J is absolutely 
continuous and tr(J') = Ui[A^_ 1; A~], where Aq < A^f ^ A^ < ... < A~ ^ Ap_ x < A~ and 
A^ are 2-periodic eigenvalues. The intervals [A^^A^n e are separated by gaps 7„ = 
(A~, A~) of lengths |7„ | ^ 0. If a gap 7„ is degenerate, i.e. |7„| = 0, then the corresponding 
segments a n , o~ n+ i merge. Recall that the Lyapunov function A(z) = an( j the 

spectrum a{J) = {A G C : A(z) G [-1, 1]}. Note that (-l) p -"A(A±) = 1 for all n G N p . 
We recall well-known facts (see e.g. [F] or |KKu| ). 

Theorem 1.1. A real polynomial F is the Lyapunov function for some scalar p-periodic 
Jacobi operator with numbers o n >0,6„6l,fi6Z iff F(z) = cz p + 0(z p_1 ) as z — > oo for 
some c > and F'(zj) = 0, (— l) p ~^F(zj) ^ 1 for all j G N p _i and for some Z\ < ... < z v -\. 

Remark. 1) Zeros of A — t for any fixed t G [—1,1] and a constant c > determine A. 

2) o~(J) = o-(J) for some Jacobi operators J, J iff A = A, where A, A are the corresponding 
Lyapunov functions. 

3) If a n — 1, b n — for all n G Z, then A = ^p(f ), where ^ is the Chebyshev polynomial. 
Before we describe the content of our paper, we briefly comment on background litera- 
ture for matrix-valued Jacobi operators. Inverse spectral theory for scalar periodic Jacobi 
operators is well understood, see a book pQ and papers |BGGK| . [K] . |KKu| . [P] . |vM| and 
references therein. The corresponding theory for periodic Jacobi matrices with matrix- valued 
coefficients is still largely a wide open field. Some new particular results were recently ob- 
tained by Gesztesy and coauthors [CG^. |CGR] . |GKM] and see references therein. Note that 
for finite Jacobi matrices with matrix- valued coefficients the complete solution of the inverse 
problem was given recently by Chelkak and Korotyaev |CK| . 

The eigenvalues of M(z) are the zeros of the equation D(r,z) = 0. This is an algebraic 
equation in r of degree 2m, where the coefficients are polynomials of z. It is well known (see 
e.g. Chapter 8, [Fo]) that the roots Tj(z),j G constitute one or several branches of one 
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or several analytic functions that have only algebraic singularities in C. Thus the number 
of eigenvalues of M(z) is a constant N e with the exception of some special values of z (see 
below the definition of a resonance). There is a finite number of such points on the plane. 
If the functions Tj(z),j G are all distinct, then N e = 2m. If some of them are identical, 
then N e < 2m and M(z) is permanently degenerate. 

If m = 1, then the Riemann surface for the multipliers has 2 sheets, but the Lyapunov 
function is entire. Similarly, in the case m ^ 2 it is more convenient for us to construct the 
Riemann surface for the Lyapunov functions given by Aj(z) = \{Tj{z) + tJ (z)), j G N m , 
which has m sheets (see the equation fll.Tj) ) . Note that the Lyapunov-Poincare Theorem 
gives <j{S) = Ujlil 2 : Aj(z) e [ — 1) !]}■ Let j^A be a number of elements of a set A. 

Theorem 1.2. There exist analytic functions A s , s = 1, ..,po ^ m on the p s - sheeted Riemann 
surface & s ,Ps ^ 1 having the following properties: 

i) There exist disjoint subsets u s C N m , s G N Po ,|Ja; s = N m such that all branches of 
A a , s G N Po have the form Aj(z) = \{Tj(z) + t^ 1 (z)), j G uj s and satisfy 

D( \ Po -I- _1 

^ ' 1 j£^s 

for any z, r G C, where the functions & a (z, v) are some polynomials with respect to u, z G C. 
If Aj = Aj for some i G u>k,j G uj s , then $fe = $ s and A& = A s . 

ii) Let some branch Aj,j G N m be real analytic on some interval Y = (a, (3) C R and 
-1 < Aj(z) < 1 /or any z eY. Then A'^z) 7^ /or eac/i 2 G K . 

All functions p,p S} s G N po ,#u; s ^ 2 ^wen 6?/ (11.81) are polynomials and 

N 

P = Y\_P S ' P^ z ) = Ft ~ A i( z )) 2 ' 2 G C > «^ ere = 1, «/ = 1- (1.8) 

w) £ae/i endpoint of a gap (A~, A+) is a periodic (or anti-periodic) eigenvalue or a real branch 
point of some Aj,j G N m; which is a zero of p (below such points are called resonances), 
v) Let Tj,j G N m be eigenvalues of a matrix A p = (a 1 a2..a p )"~ 1 . Then following asymptotics 
hold 

A j {z) = j{7f + 0{z-*)), T j (z) = 2?Tf + 0{z-±)), j G N m , (1.9) 
p s ( z ) = ^l(l^|-l)( Cs + 0(^)), c s = 2^K 1 -I- S |) JJ (rO-rO) 2 , s G N Po (1.10) 

as \z\ — > 00 Moreover, if rj 1 7^ r° /or a/i j, fc G cj s , j 7^ fc, t/zen c s 7^ and p s {z) = 
z | w .|(| w .|-i)( Cs + 0{z- 1 )) as \z\ -> 00. 

Definition. number z is a resonance of J , if z is a zero of p given by (II .81) . 



3 



Example. Let a n = diag{ai jn , a 2 ,„, a m ,„, } and b n = diag{&i jn , b 2 - n , — , b m>n , } for all n G Z 
and for some a^ n > 0, bj >n G R, (j, n) G N m x Z. Then the operator J7" = @™Jj, where J} is 
a scalar Jacobi operator, acting in £ 2 (Z) and is given by 

(JjV)n = aj,nVn+\ + aj,n-iy n -i + bj, n y n , n G Z, y = (y n ) ne z G £ 2 (Z). 

In this case each Aj is the standard Lyapunov function for and the properties of Aj is 
well known, see Theorem 11.11 and [Tj . Thus p = 1, since cjj = {j}. ■ 
Below we show the following identity (see (13.21) ) 

mp 

D(*,t) =cr m Y[(z- Xj(r)), z, r G C, r ^ 0, c = (-l) m det A p , (1.11) 

3=1 

where A 3 -(r) are zeros of the polynomial D(z,t) = for fixed r. Define the set A(r) = 
{A n (r), n G N mp }, r G C. Then we obtain <r( J) = [jf^ a jt where ^ = A^S 1 ), S 1 = {r G C : 
|r| = 1} is the bounded close set. Ler = {0, ...,m} for m ^ 0. We describe minimum 
spectral data which determine all Lyapunov functions Aj,j G N m for some Jacobi operator. 

Theorem 1.3. Let x fc G satisfy cosx n ^ cosxj for all j ^ n. 

i) LetA = A(e l *°) and let Aj C A(e^ J ), #Aj = (m—j)p+l,j G N m for some Jacobi operator 
J. Then the spectral data Aj,j G N^, determine the polynomial D(-,-), all Lyapunov 
functions Aj,j G N m and £ne spectrum cr(J') (counted according to its multiplicity). 

ii) Let Ai C Ai and let #Ai = t^Ax — 1. TTien inere exist infinitely many Jacobi operators, 
having the same spectral data Aq, Ai, Aj, j = 2, ...m, but different determinants. 

Example. Note that if e lk>tj = 1 for some k G N, then X n (e l><j ),n G N mp are k-periodic 
eigenvalues. Consider the case m = 2. Let x = 0, xi = ir. Then A is a set of all periodic 
eigenvalues and Ai^is a set of some anti-periodic eigenvalues such that #A = 2p, #Ai = p+1. 
Let the set A 2 = {A}, where A G o~{J) is a some four-periodic eigenvalue, i.e., x 2 = f . Due 
to Theorem 11.31 the spectral data A ,Ai,A 2 determine the polynomial det(^f p (z) — rl m ), 
the spectrum cr(j7") and the Lyapunov functions Ai, A 2 . ■ 

Theorem 1.4. i) Let p ^ 3. Tnen eac/t swm X]j=i Aj(r), s G N p _i does noi depend on 
t G C. In particular, the following identities hold 

mp p mp p 

X> n (T) = ^Tr6 n , E^M = ^Tr(6 2 +2a n aI), a// r G C. (1.12) 

n=l n=l n=l n=l 

nj The following estimate is fulfilled: 

pm 

A^(r) ^ 2pm(det A p )^, (1.13) 

n=l 

2 

where the identity holds true iff b n = 0, a n a„ = (detA p )p m I m for all n G N p . 
in) Let numbers Xj G M,j G N m satisfy cosx n 7^ cos Xj /or a// j 7^ n. Then the eigenvalues 
A nm+fc (e i ^) = 2 cos i(xj + 2im) for all (j, n, k) eN° m x N°_ x x N p (ff K = 0, a n aj = L m for 
all n G N pm and EKUi a n = I m . 
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We describe a priori estimates 

Proposition 1.5. Let Ht/Hoo — max {l a n(j) k)\, \b n (j,k)\} } wherethe matricesa n = {a n (j, k)}, 
b n — {b n (j, k)}. Then the following estimates hold 

\\JWoo ^ max{|A+|, |A^|} = \\J\\ s= (4m - l)||^T||oo, (1-14) 

lUlloo + ' AjV \ A ° +l < < (4m - \)\\J\U if £>&,- = 0. (1-15) 

j=i 

Remark. A priori estimates for scalar Jacobi operators were obtained in |BGGK| . [K] . 
[KKr] . 

The plan of our paper is as follows. In Sect. 2 we prove Theorem II .21 In Sect. 3 we 
prove Theorem 11.31 11.41 and Proposition 11.51 In Section 4 we consider examples for the case 
m = p = 2, where we construct the complex and real resonances. In Sect. 5 we shortly recall 
the well known results about the properties of point spectrum and the absence of singular 
continuous spectrum. In the proof we use [CKJ, [BBKJ. 



2 The Lyapunov functions 

We recall well-known results. For any solution y = (y„) ne z of the equation a n y n+ \ + b n y n + 
o7 n _ x y n ^x = zy n we define 

= ( £ ) ' ^ X = ( a-A^- K) ) ■ " £ Z ' ^ 

Then f n satisfies f n = T n f n -\. Thus the matrix- valued function = ( f n r " 

satisfies M n = T n ^ n _i, ^ = I 2m - This gives the monodromy matrix *4t n = Y]^ = i7j, 
n ^ 1, where rij=i^i = X n ...X\ for matrices Xj. We rewrite T n in the form 

r«=('z £ ) ( _] m , i\ ) ( °r I ) - («) 



,T m T D D / Im \ T m -1 t-0 ' -^n 



where 

P n = 0>n © ^m> -Rn — -?n ( J* J = a n © ^ = ( _j ^ J ■ (2-3) 

Then using (12.21) . we get 

n n n 

^n = H I) I = P" 1 (II P = Pn'MnPo, M n = J[ R/lf '. (2.4) 

j=l j=l j=l 

Matrices Rj,7j, j ^ 1 satisfy 



RJJR^J, (T°) T JT° = J, J=(_° Im ^) 



(2.5) 



M n = P n M n Pt = ( i^Hll ^ ), Ml={ a °T n+1 ) (2.7) 



which yields 

.\/„' J.\l„ = J, M" 1 = -JMjj, n G Z. (2.6) 

Using ( 12.41) we obtain 

dn+ii^)- 1 <f n+1 J ' ±v±n ~ \ <p T a n (fl +1 
and 

M" 1 = -JMjj = ( ~_f^ n ) . (2.8) 

Due to .\l„.\l„ ' = /2m and M~ x M n = I 2m we deduce that 

<£0„(4)-V£fi " <£Pn(4 = U i?„(aJ)-Vn = ^(4) _1 C (2-9) 

4<WrH-i - 4+i a nV?n = oq, (p~Z +1 aZp n = ^aj^n+i, 4+i a T^n = $Ja^ n+1 . (2.10) 
Recall the asymptotics of fundamental solutions 

<p n+1 = z n A n + OiX^ 1 ), tf n+1 = O^- 1 ), n ^ 1, A n = (oxaa.-On)- 1 (2.11) 

as z — > oo. Substituting last asymptotics into we obtain 

^ P W = ^(o M+O^" 1 ). (2.12) 
It is well known that the determinant D(z,t) satisfies 

2m 2m 

D(z, r) = det(M(z) - rl) = r 2m + £ ^(,(2) = JJ(r - 73 (*)), a, r G C, (2.13) 

3=1 j=i 



£o = l, 6 = -Ti, 6 = ~(T 2 + T^ 1 ), X,.- (2-14) 

^ s 

see p. 331-333 [RS]. The identity (11.60 gives £j = (,2m- j for all j G N m . Recall that the 
Chebyshev polynomials <%,, n ^ 1 satisfy: 

f-1 

W = T —^~ = 2™- 1 £c^-« = (-!)'" (n-V)!^ 2 "" 2j ~ m ' (2 - 15) 

i/ = 2^f^, see [AS]. Then ( 12131) and the identity = yield 

^ = + 6^1^> + ... + 2-6. = ± 

V ^ 3=0 
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1 



TiM n (z) 



^ T - + T- 



2 v ' ^ 2 

1 

The substitution of (12 . 151) into the equality (|2. 161) gives 

D(r,z 



3=1 



z e 



(2.17) 



(2r)' 



2 1/ 



(2.18) 
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! Cm-l,oCl — 7T) 02 — C TO| i + C m _2,0^2, 03 — Cm-l,l6 + c m-3,0^3, •••) (2.19) 
02n — Cm,n + Cm-2,n-l£,2 + C m _4 jn _2^4 + ••• + Cm-2n,o£,2n, (2.20) 
02n+l = Cm-l.raCL + C m _3 ;n _i^3 + 2^5 + ••• + C m _2n-l,0^2n+l- (2-21) 

Let a n = I m ,b n = for all n G Z and denote the corresponding Jacobi operator by J° . If 
/}° is the corresponding determinant, then we obtain D°(r,z) = (r 2 + 1 — 2r^(|)) m and 

" ~ iSjS* = (" - ^(D) m = E"^(-^(f))i^, where = 



Proof of Theorem 11.21 i) We have the identity $( 



Em 




see 



( 2r )m A^O Yj\ 

(12. 181) . where the polynomials 0j are given by ( 12.1 9h - f T2~T2~TT) . Ai(z), .., A m (z) are the roots of 
$(z, i/) = for fixed z G C. Then the statement i) follows from the well-known about the 
zeros of a polynomial ^™ v m ~^uij (z) , where iWj is a polynomials in 2, see |Fo| . 

ii) We have A'- (2) = |(1 — r _2 (2;))r / (2;) 7^ 0, z G F, since by the Lyapunov-Poincare 
Theorem, r^z) 7^ for all zGF. 

iii) Recall that the resultant for polynomials / = r n + a\T n ~ x + .. + «„, g = &o rS + 
l3\T s ~ l + .. + p s is given by 



i2(/,»)=det 
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(2.22) 



The discriminant of the polynomial / with zeros Ti, .., r n is given by 

Di S /=n^-^) 2 = 



n(n — 1) 

-l)^R(fJ'). 



i<j 



_ iV(JV'-l) 

Thus we have Dis^(A,r) =Ui< s ,i,se4 A i( X ) ~ A s(X)) 2 = (-1)— fl^,^-) is polyno- 
mial, since $j(A, r) is the polynomial. Then the function p = Y\i° Dis $j is polynomial. 

iv) Each gap has the form 7„ = (A~,A+) = n^ =1 7 n j, where 7 nj - = (A~j,A+j) C K is 
some finite interval such that Aj(z) £ [— 1, 1] for all 2 G 7 n j. Note that Aj(A^j) = ± or 
Aj(A^j) = ±1 or A„ • is the branch point of Aj(z), otherwise we have a contradiction. 
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v) Using (1212]) we define F by 

F(t) = t p ^ p = Fo + tF^t), F = O m ®A p , t = z-\ 

where F\ is a some matrix polynomial and F\{t) = const +0(t) as t — > 0. Let Tj(t),j G N2 m 
be the eigenvalues of F(t). Recall that t®, j G N m are eigenvalues of A p . Then, Tj(t) are 
given by convergent Puiseux series (the Taylor series in r 1 /"^ ) 

°° k ^ 

Tj(t) = T® + ^ OLj^V n i , \t\ < r, mj = rn, mj G N, j G tu S) s G N fe 
fc=i 1 

for some r > 0, see p=4, [RS], where W\ = N mi , w 2 = N mi+m2 \ u>i, u; fc = N m \ N m _ mfc 
and rj 1 = r° for all j, j' G w s . Here Tj(t),j G u> s are the branches of one or more multivalued 
analytic functions with at worst algebraic branch points at r = and satisfy 

Tj{z) = Z P Tj(t) = Z P (t° + 0(z~^)) = T°Z P + 0(z P ~^), J G W s . 

If each eigenvalue rj 3 is simple, then each Tj(t),j G N m is analytic for small t, i.e., m s = 1. 
Thus we obtain Tj(z) = t°z p + 0{z p ~ l ),j G N m , which yields 2Aj = Tj(z) + ^ = t°z p + 

0(z p ~ 1 ) and using (II .81) . we get asymptotics for p s . 

3 Inverse Problems 

In order to prove Theorem 11.31 we need the Lemma 3.1-3.3. 

Lemma 3.1. Let W = {W r j}j r=Q be a s x s matrix with components W r j = cosjx r for 
some x r G C, s ^ 1, where cos x r 7^ cos K r > for all r 7^ r' . Then det W 7^ 0. 

Proof. Using cos jz = ^(cosz), where the Chebyshev polynomials ^ are given by (j2.15l) . 
we obtain W r j = cos jx r = ^ (cos x r ). Let Wy = for some vector y = (y n )o £ C s+1 , y 7^ 0. 
Then we deduce that a polynomial P(x) = ^2j =0 yj^j( x )i P 7^ has distinct zeros cos x^, j G 
N°, which gives a contradiction, since degP ^ deg S? s = s. Thus W is invertible. ■ 

Lemma 3.2. Let c = (— l) m detv4 p . The polynomials £j(z),j G N m , given by (j2.13l) . satisfy 
Zj{z) = bm-Az) = 0{z p i), U(z) = cz pm + O^ 1 ) as z - 00, , (3.1) 

mp , . m— 1 



n=l j=l 

Proof. Using ( ll.6p . (12.1 3H and Viete's formulas, we get 

j 

Uz) = bm-Az) = il^W' (3-3) 



Asymptotics ( 11.91) give Tj(z) = TjZ p + o(z p ), j G N m and the Lyapunov-Poincare Theorem 
yields T m+ j(z) = ^-^y = 0(z~ p ) as z — > oo, since each ^ 0. Then, using (13.31) . we obtain 

m 

Zj(z) =0(z P1 ), U{z) = (-l) m Y[r°z pm + o{z pm ) = (-l) m det A p z pm + 0{z pm - 1 ). 

3=1 

Identity follows from last identity and (Oil . (12T3D . ■ 
We need the following Lemma about polynomials. 

Lemma 3.3. Let (z n ) k G C k and (cj)q G C s+1 /or some integers s, k ^ 1. TTien £/iere eziste 
a unique polynomial r(z) satisfying for some polynomial g the following relations 

s k 

r(z) = z k ^2cjZ j + 0(z k ~ l ) &sz->oo and r(z) = h(z)g(z), h(z) = JJ(z-z n ), (3.4) 

j=0 n=l 

Proof. Introduce the linear space of polynomials V s = {g : degg ^ s}, dimP s = s + 1. 
Note that g E V s . Define the linear operator A : V s — > C s+1 by 

1 rf n / \ 

(^)j = ^y— (M*)#(*)J | 2=0 , n = A; + j - 1, j G M s+1 , # G 7> a . 

We rewrite Ag in the form (Ag)j = Cj_i, j G N s+ i, where h(z)g(z) = z k ^2j=o CjZ-' + 0(z k ~ 1 ) 
as z — > oo. We deduce that ^4 : "P s — > C s+1 is an isomorphism, since dimP s = dimC s+1 = 
s + 1 and if Ag = 0, then g = 0. If g = ^4 _1 (c n )Q and r(z) = h(z)g(z), then r(z) = 
z k ^2 S j =0 CjZ^ + 0(z k ^ 1 ) as z — > oo. Moreover, the polynomial r is unique, since A is an 
isomorphism. ■ 

Proof of Theorem II. 3L Recall that c = (— l) m det A p . Define the polynomials (o(z) = c -1 , 
Q(z) = c _1 ^j(z), j G N m . Lemma f3T2l yields degCj ^ pj and then Cj(z) = Y^=oCj,nZ n fo r 
some Cj,n G K. Introduce the sets -?T S = {p(m— s— 1)+1, ...,p(m— s)}, s G N^_ 1; K m = {0}, 
which satisfy 

j s= j j\ s = i^ p (- m _^, 
Using (13.21) and (13 .5p and setting hj = T~i + r J , we obtain 



UZiKs = K(m -4), K r nK, = Qjor all r ^ s. (3.5) 



D(z r) pm 

l( Z ' T ) = CT L = 11^ ~ A ™( T )) = Yl h m-jQ( Z ) = Yl h 3<m-j{z) (3.6) 



n=l j=0 j=0 



m p( m ~j) mm mm 

fy? Cm-j,nZ = /lj C m _j )rl 2 = ZljC 

j=0 n=0 j=0 s=j n£K s j=0 s=j n£K a 



m—j,n,Z 



m s m s m pm 

s=0 j=0 n&K a s=0 n&K a j=0 s=0 n€K a n=0 
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where 

s 

Vn(r) = h j(m-j,n, neK s , se N° m . (3.7) 

3=0 

Substituting r = e l><r ,r G N° into the identity ( 13.71) we get 

?j n = 2W s y n , rjn = {Vn{e i><r )) s r=0 , y n = (Cm-j,n) s j=0 , W s = (cosjx r )* i=0 , (3.8) 

for all (n,s) G K s x N^. Due to Lemma [37TI each matrix W s ,s G N° m is invertible. Thus 
0, (132!) give 

r/n(r) =< (r j + r^)^, (2Wg"^ n >, (n, s) G K s x (3.9) 

where < •, • > is a scalar product in M s . Then the functions T) n (-),n G are determined by 
Vn — (Vn( el:>er ))r=o- We w ^ use this fact to determine all functions T) n (-),n G N m . In order 
to describe recovering we need the following simple fact. 

Lemma 3.4. Let < k < m - 1. T/jen numbers r] n (e l: ^), (n,j) G N° m x N° ; £/ie /unc- 
tions n G U?? =0 i£j, t/ze set A k+1 and the numbers x,, j G determine the numbers 
Vn(e iitj ), (nj) ^ N° m x Ng +l3 and tne functions rj n; n G U^if,-. 

Proof. Let S k+1 = (^(e^OWu^ = fa«(e^ + 0)S,( m - fc -i )+ i- Using (ED, we see 
that the elements of A fc+1 are some zeros of the polynomial q(z, e li<h+1 ) and recall that 
#A fc+1 = p(m — k — 1) + 1. Also, using (13.61) . we see that the components of vector Sk+i 
are the coefficients of the polynomial q(z, e l>Ck+1 ). Then, due to Lemma 1331 A fc+1 and S k +i 
uniquely determine the polynomial q(z, e l>Ck+1 ) and, using (13.61) . we determine (r] n (e l:><k+1 ))^ . 
Substituting obtained values (T)n( eti<:i ))neK h+1 , j G N° +1 into the ( 13.91) . we determine r} n {r), 
n G K k+1 . * 

Now we describe the recovering procedure. Recall that q(z,r) = D<sZ £) . Substituting 
the elements of the set A into the (13.61) . we determine q(z,e l>co ) and (?7ri(e iJ<0 ))^™ . Then 
( 13.91) determines all functions r) n (r),n G K . Step by step, using Lemma [331 for k G N^_ l5 
we determine r) n (r), n G {J-L Q Kj = Np m . Then (13. 6f) gives q(z,r). Using (12.131) we get 
cq(z,r) = D(z,t) = r 2m + 0(r 2m_1 ) as r — > oo, then the polynomial q(z,r) uniquely 
determine constant c and D(z,r). 

ii) Assume that a n ,b n ,n G Z are diagonal, and each a n > 0,n G Z. Then J = 
where are scalar Jacobi operators. Let Aj be the Lyapunov functions and let Dj be the 
corresponding determinant for the operators G N m . Then we deduce that 

m 

D{z,t) = Y[D v Dj(z,r) = t{t + t- 1 -2A j {z)). (3.10) 

3=1 

By Theorem 11.11 the Lyapunov function A m (z) satisfies A' m {z®) = and (— l) p ~ J A m (,2°) = 
|A m (z°)| ^ 1, j G N p _i for some z\ < ... < z®. Then the polynomial q s (z) = s(A m (z) — 
cos x ) + cos xo, s ^ 1 satisfies 

(-l) p -^ s (^) = s(|A m (z°)| + (-l)^(-l + i) cosx ) ^ 1, ( ?s )'(z°) = 0, j G N^. 
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Using A m = c m z p + 0(z p ~ 1 ), c m > as z — > oo, we get q s (z) = CmSz^ 1 + 0(z p_1 ), c m > 0. 
Due to Theorem 11.11 q s is a Lyapunov function for some scalar p-periodic Jacobi operator 
J^. Consider the p-periodic m-dimensional Jacobi operator J s = © Jj. Then the 
determinant for J s is given by 

m— 1 

D s (z, r) = r m (r + r' 1 - 2q s {z)) JJ (r + r' 1 - 2A i (z)). (3.11) 

3=1 

Using (I3TTUL pTTTil . we see that A (J r ) = Ao(.7 s ), where Aj(J),j = 0, 1, 2.. is the set Aj for 
J . Each Lyapunov function Aj 1 (z),ji G N m ._i satisfies 

p 

A h (z) - cosxj = cjj JJO - z n+pil j), Cy-j > 0, j G N m , for some z n+pjlJ G R. 

n=l 

Then using (13.1 Oft . (13.111) . we deduce that z n+p j,n G N p ( m _i) are zeros of each polynomial 
D s {z,e i3 ^),j 6N mi O 1, where D x (z,t) = D(z,r). We take Aj(J s ) = Aj(J) = {z p+n ,n G 
N p („_ j)+ i}, j = 2, ..,m and Ai{J s ) = A\{J) = {z n+p>1 ,n G N p ( m _i)}, where A\{J) is the 
set A x for J . But we have D s ^ D\ = D for any s > 1, since q± = A m ^ q s . ■ 
Proof of Theorem Q i) Let g n (r) = (-l)"£ ls£il< .. <Jn<mp IEU A^r). Lemma 
implies = 0(z PJ ). Then using (13.21) for fixed r 7^ ri = 1, r 7^ 0, we obtain 

cf^z^Mn) - ffi (r)) = cj](z - A^n)) - cfl(z - Ai(r)) = 

m— 1 

= W + - r m - r~ m ) + ^(rf-" + r"- m - r m - n - r n " m )e„(^) = 0(;z (m ~ 1)p ), 

n=l 

which yields gj{r\) = <?j(t) for all j G N p _i, since c 7^ 0. The Newton formulas give 

n—l mp 

ui = -g h u n = -ng n - Efl^n-j, where M n (r) = A"(r), 

3=1 3=1 

which yields that m„(t) = « n (ri) = u n {l) for all (n, t) G N p _i x C. Due to (15.41) . we get 

mp mp N 

£aJ(t) = TrL(rr,0 1, E A ^ = TrL ( r ) = E Tr ^ 
3=1 3=1 3=1 

mp p p 

£ A|(r) = TrL(r) 2 = TrL(r)L*(r) = 2^ INI' + E INI*- ( 3 - 12 ) 

3=1 3=1 3=1 

ii) Let £j,ns J ^ N m be eigenvalues of the matrix G Np. Then ( 13.121 ) gives 

pm p pm m _j_ 

E A «( r ) ^ E Tr2a « a « 2 EE > ^ 2 P m (n^) Pm = 2pm(det^)^. (3.13) 

n=l n=l n=l j=l j,n 
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Recall that the identity Y7n=i Sj=i £ j,n = P m [Ylj n £ j,nj P holds true iff each e^™ = det A 2 p , 
(j,n) G N m x N p . Then the statement ii) holds true. 

iii) Sufficiency. Let a = (a°) ne ^, b° = (6°) nG z, where a° = I m , 6° = for all n G Z. Let 
the Jacobi operator J" = J7"(a, b°), where each matrix a n , n G N p is unitary and a n = ^m- 
Define matrices c n by Ci = J m , c n+ i = a^c n , n G Z. Thus c n+p = c n for all n G Z and c n are 
unitary matrices, since a n are unitary. Then £/ = diag(c n )„ 6 z is unitary operator in £ 2 (Z) m 
and £7,7(0, fr )^ 1 = J(a°, b°) = J°, which yields 



in 



det(M p - rl 2m ) = r m JJ(r + r- 1 - 2& p (z/2)) = det(M p ° - rl 2m ), 

5=1 

where M p , M° are corresponding determinants and 2? p is the Chebyshev polynomial, see 
(12.151) . Using ni.lll) and ^£(2) = cosQoarccosz), we obtain \ nm+ k(e l>Cj ) = 2cos-(xj + 27m) 
for all (j, ra, fc) G x N°_ x x N p . Also, using i), we obtain 

pm p 

53 An(e^) 2 = E Tr ((^) 2 + 2 « T ) = J G N - ( 3 - 14 ) 

n=l n=l 

Necessity. Let A nm+fc (e^) = 2 cos ~(x,- + 2?m) for all (j, n, k) eN° rn x x N p . Then 
Theorem 1 1 . 3L ii uniquely determines 

m 

r' m D(z, r) = J](r + r- 1 - 2& p (z/2)) (3.15) 
i=i 

(see above Sufficiency). Identities ( 11. lip give detA p = 1. Using ( 13. 14h and Corollary ll.41 ii. 
we obtain b n = 0, a n a^ = I m . Due to ( 11 .7L d 1 - 9 [) and (13.151) . we obtain 

A, (z) = l{ Tj (z) + rr\ z )) = T Jf + O(z^), Aj(z) = %{z/2) = | + O^ 1 ), 

then rj 1 = 1, j G N m , i.e., all eigenvalues of matrix A p are equal 1 and all eigenvalues of 
matrix A' 1 = a\...a p are equal 1, then A p = J m , since a n (a n a^ = I m ) and are unitary. ■ 
Proof of Proposition [TT5l For the self-adjoint operator J we obtain | |J7"|| = max{|AQ |, |A |}. 
Moreover, the operator J has the form J = "Y^^m-i^h wnere •Jj is some diagonal matrix. 
Note that Jj,J s ,j 7^ s have different diagonals. Using the identity \\Jj\\ = WJ'jWoo and the 
simple estimate ||i7||oq ?5 ||i7|| we obtain 

2m-l 



UlU < ||t7|| ^ Yl (4m -1)11.71 



-2m+l 



which yields (fTT4lh If E? Tr& i = °> then £Aj-(t) = 0, which yields A^ < ^ A and 
A — Aq = |Aq I + |A^ |. Thus all these stimates gives (11.15p . ■ 
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4 Examples for the case p = 2, m = 2 

We consider the case p = 2, m = 2, where real 2x2 matrices a n , b n satisfy 

/ 1 02n+l \ 7 _ 7 / «2n /?2n \ _ /3 _ /3 

«n+2 — a « — I n 1 ' °™+ 2 ~~ °™ ~~ I a „ ' a n+4 — &n, Pn+4 — /V 



.0 1 ; ' "™ V &n «2n+l . 

The determinant D, the function p and the Lyapunov functions are given by 

£>(-, r) = r 4 - Tir 3 + Tr 2 - 7\r + 1 = (r 2 - 2Axr + l)(r 2 - 2A 2 r + 1) = (2r) 2 $(-, i/), 

T _T[-T 1 _ -T 2 + 2T 2 + 8 _ T 1 -^p _ 7\ + y/4p 
_ 2 ' P " 4 ' Al " 4 ' A2 " 4 ' 

T 1 = 2z 2 + T 11 z + T 10 , T 2 = z 4 + T 23 z 3 + .. + T 20 , $(•, v) = v 2 - ^-v + j - ^ 

where 

3 

T n = - «n + (A) + + A(A) + (h), 



3 

Tio = -A + 2p p 2 + 2p 1 p 3 + ]Jp n + aoa 2 + a 1 a3-Po(a2P3 + Oi3Pi)-P2(aoPi + ai/3 3 ), (4.1) 

o 

3 3 

t 23 = - ^2 an > T<22 = ^2 a i ak ~ ^2 & ~ 4 ' 

0<j<fc<3 

?2i = (a 2 + a 3 )(3l + (a + a 3 )0l + (a + «i)/3 2 + («i + a-zjPt 
— 2/3 (/5i + A$) - 2/5 2 (/3i + /3 3 ) + («i + e* 3 )(2 - a a 2 ) + («o + a 2 )(2 - aia 3 ), 

3 

T 20 = 6 + /3 2 /5 2 + /3 2 /5 3 2 - 4/3^3 - 4(3 (3 2 + JJ a n - a (ai/3| + a 3 /3 2 ) - a 2 {a^ 2 3 + a 3 /3 2 ) 

o 

-2a a 2 - 2 ai a 3 + +2(3 (a 2 (3 3 + a^) + 2 / 3 2 (a 0/ 3 1 + Ms)- (4-2) 
Below we consider some specific cases. 

1. All a n , b n are diagonal matrices. In this case all (5 n — 0, the function p = and 

1/ Q + a 2 \ 2 l/ ct -ct 2 \ 2 1 / Qi + Q 3 \ 2 l/ ai-a 3 \ 2 

1 ~ 2 V 2 / 2 V 2 / ' 2 ~2V 2 J 2 V 2 J 

2. The constant coefficients. Let a n = 0, (3 n = f3 for all n G Z. Then 

Ti = 2^ 2 + 4(3 2 z + [3 4 + 4[3 2 - 4, p = /3 2 (2z + (3 2 ) 2 (4z + [3 2 + 8). 

2 

The function p has 3 zeros z±, z 2 , z 3 . If (3 — 0, then p = 0, Ai = A 2 = y . 
If /3 ^ 0, iv 7 ^, then p has the zeros z x = z 2 = z 3 = 
If /? = ±\/8, then p has the zeros z\ = z 2 = z 3 — 4. 



- 1. 
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We have the following equations for periodic and anti-periodic eigenvalues 

4(Ax - 1)(A 2 -l) = (z + 2) 2 ((z - 2) 2 - A(3 2 ) = 0, 4(A X + 1)(A 2 + 1) = (z 2 - 2(3 2 ) 2 = 0. 

Hence for each (3 ^ there exist two double anti-periodic eigenv alues -/V2,/V2; one 
double periodic eigenvalue —2 and two simple periodic eigenvalues 2(1 — b), 2(1 + b). If 
j3 = ±2, then there exists the periodic eigenvalue —2 of multiplicity 3. 
3. Example (a , «i, a 2 , a 3 ) = (1,0,-1,0); (fo, Pi, fa, (h) = (*, 0,0,0), t 6R \ {0}. In this 
case we have 

T 1 = 2z 2 -5, T 2 = z 4 - (t 2 + 5)z 2 - t 2 z + 8, Ap t = At 2 z 2 + At 2 z + 1. 
The Lyapunov functions and resonances rf (zeros of p) are given by 



2z 2 - 5 - jAp t 



A 2)t = 



2z 2 - 5 + v/4p^ 



-l ± v / T 3 ^ 



Each resonance , £ G (0, 1) is complex. 

Let t — 1. Then the Lyapunov functions are given by A 1;1 = 



,A 2 ,i 



™^ and 



the spectral bands a\ = {z G R : Aj,t( z ) £ [ — 1) 1]}? J — 1? 2 have the forms 
(7 1 1 = [Ai,Aj]U[Aj,Ai] = 



v^T-1 v 7 ^- 1] rv / 5 + 1 v^+l 



^[A^JU^A 2 



2 2 
v / 17+ 1 



-1 



U 



0, 



2 ' 2 
v / 17+ ll 



where A^j = 1,2, k = 1,4 are periodic eigenvalues, Aj,,j = 1,2, k — 2,3 are anti-periodic 
eigenvalues and the resonances rf = — \ lies on the gap (—1,0). 

Let t — 1 + e for some small £ > 0. Then the spectral bands have the forms 

^ = [A( t ,Ai t ]U[Ai )t ,Al t ], Aj t = Aj + o(l), yeN 4 xN 2 ,r t ^^ as t | 1. 

Note that the "resonance" gap j t = (r t ~, rf) 7^ 0, since rf < rf and cr(J) D 77 = as t I 1. 
Thus the interval 7$ C (A| t , A 2 1 ), where (A 2 , t , A 2 1 ) is a gap in the spectrum of J . This gives 
that "resonance gap" arises in the gap and the end points of the spectral bands are periodic 
or anti-periodic eigenvalues. 

If t = 1, then the resonances rf = rf are not the branch points. If t > 1, then the 
resonances r\ < are the real branch points. If t G (0, 1), then the resonances rf G C± 
are the complex branch points. 

Example 4. Let (a , a u a 2 , a 3 ) = (0,1,0,1); (fy, fa, fo, As) = (t, 0,0,0), t G R. Then we 
obtain 

T x = 2z 2 -2z-3, T 2 = z 4 - 2z 3 - (3 + t 2 )z 2 + (4 + t 2 )z + 4, 
Ap t = (2z - I) 2 + At 2 (z 2 -z), rf = \ T G R. 
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Here r t are the zeros (resonances) of pt and if t 7^ 0, then r t < . Thus we have the 
resonance gap ^ t = {r\ ,rf) C (0, l),t 7^ 0. The Lyapunov functions are given by 

Ai, t = l -{2z 2 -2z-3- y/4p t (zj), A 2)t = -^Lz 2 - 2z - 3 + ^(2)). 

Let t — 0. Then the Lyapunov functions have the forms Ai = 2 ~ f^— K A 2j0 = ^-f^ 
and A lj0 (|) = A 2;0 (|) = -|. The spectral bands of = {zG»: A j:t (z) G [-1,1]}, j = 1,2 
have the forms 

(To 1 = [-1,3], <x 2 = [-2,2]. 

The end points of spectral bands <Tq, ct 2 are periodic eigenvalues. 
If t G (0, e) for some small e, then we obtain 

^ = [Ai(t),r t -]U[r+Aj(0], r~ < r+ rf = ^il, 

Aj jt = -l + o(l), A^ = 3 + o(l), A? jt = -2 + o(l), A| t = 2 + o(l) as i - 0. 

This yields cr(j7") = cr* U of = [A 2 ;t , r t ~] U [r+, A4 J, where (r t ~,r t + ) 7^ is the resonance gap, 
since r^T < r+, and A 2 1 , A 2 1 are periodic eigenvalues. 



5 Appendix 



Let Jff = L 2 (j0, 2?r), g, ft) = /® 27r) ftg be a constant fiber direct integral, where ft = C pm . 

It is the Hilbert space of square integrable ft-valued functions. Let ^(Z) be the set of all 
compactly supported functions (/ n ) ne z £ / 2| ""' 



(5.1) 



Lemma 5.1. T/ie operator U : / 2 (Z) m — > J^ 7 , given % 

00 

= E e ""V, + n P , (*, j) G [0,2tt] x N p , 



«5 we// defined for J^(Z) anc? uniquely extendable to a unitary operator. Moreover, 



UJU 



J[o. 



2tt) 



L(t) 



/ 61 ai 

&2 a 2 

2 T &3 



6, 



y ra p 



p-i 








v(J) = v a c(J) = IJ a(L(e fa ))= |J u£™{A n (0}. 

a;e[0,27r) xe[0,27r) 



(5.2) 



(5.3) 
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Proof. We use standard arguments from |RS| . For / £ J5^(Z) the sum ( 15 . 1 h is clearly 
convergent. For such functions / we compute 



\\uf\\ 



/ \\(uf)(x)\\ 2 ^= / E( E e "^w E ' 

^ j'=l n=— oo s=— oo 

P OO OO „27T J P OO 

E(E EW J+ »/ = E E ii/: 

7=1 n=— oo s=— oo ^ 7=1 n=— oo 



/ ;+np|| 



j'=l n=— oo 

Then £/ is well defined and has a unique extension to an isometry. To see that U is onto 
we compute U*. We define 

/" 27r r/r 

(U*g) j+pn = J e"^)- eC m , (j»eN p xeZ. 
A direct computation shows that this is indeed the formula for the adjoint of U. Moreover, 

ii^ii 2 - t^i E E ik^wii 2 = E E | / ' 

j'=l n=— oo jr'=l n=— oo 

= E/ ii*(«)ii'| = / Eii«.u)ii'| = ii#, 

where we have used the Parseval relation for the Fourier series. 

We verify (I5.2p . Using a n+p = a n , b n+p = b n we obtain for / £ J^(Z) 

oo oo 

{U <H f)j(%) ^ ^ ^ (J f)j+pn ^ ^ e iflj+pn—lfj+pn—l ~\~ bj+pnfj+pn 4" ®>j+pnfj+p i > + 1 

n=— oo 



~~ a j-i e mx fj+pn-i + fej e mx fj +pn + +Oj e inx fj +pn+ i 

n=— oo n=— oo n=— oo 

= a,- ,(/•/),(./ - 1) + bj(Uf) x U) + «.,(('. f)A.,). J £ N p , 
where we define (Uf) x (0) h (Uf) x (p + 1) (since (Uf))x(j) is defined for j £ N p ) by 

oo oo 

(f//),(0) = £ e- m Vpn = e"*- £ e-^-^^p+p^i) = e~ ia: (£/J%(p), 

n=— oo n=— oo 

oo oo 

(Uf) x (p+1)= e~ mx f p+1+pn = e™ e- iin+1)x f 1+p{n+1) =e"(Uf) x (l), 

n=—oo n=—oo 

which yields (UJf)(x) = L(e ix )(U f)(x). 

The eigenvalues A„(e i:r ),n £ N mp are piecewise real analytic functions on [0, 2ir] , see 
[RS]. This yields a sc {J) = 0. Moreover, the standard simple arguments from |DS| yield that 
the operator J has hot have eigenvalues. Then standard arguments from |RS| yield (15.31) . ■ 
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Lemma 5.2. The following identity hold true 

, v pm 

det{L{r)-zI pm ) = ^^ = {-iy m \[{z-\ n {T)\ all r^O,z. (5.4) 

CT 

n=l 

Proof. Let L(r)f = A(r)/ for some eigenvalue A(r) and some eigenfunction / = (/ n )n=i e 
C pm . If /o = T _1 / p , then the definition of the matrix ^ v gives ^f p (\(r))(f , fi) T = 
T (fo, fi) T - Thus r is a multiplier of J% p (\(t)) and A(r) is a zero of D(-, r). 

Firstly, let all eigenvalues X n {r),n G N mp of L(r) be distinct for some r G C. Then 
A„(r),n G N mp are zeros of D(-,t), which yields (I5.4p for all r, z G C, since the orders of 
det(L(r) — zl mp ) and D(z,r) are coincide. 

Secondly, consider the general case. Let r n = n,n G N mp and let r n+pm = r n for all 
n G Z. Define the operator Jacobi Jt = J + t diag(r n ) ng ^, t G IR and the corresponding 
matrix L t (r) = L{r) + 1 diag(r n )^™ 1 (given by (15.21) ) and let D t (z,r) be the corresponding 
determinant. Then all eigenvalues A t n (r),n G N mp of L t (r) are distinct for some r G C and 
t — > oo, since A tjn (r) = tr n + o(t) as t — > oo. Thus we obtain 

D ( ) pm 
det(L t (r) - ^ pm ) = = (-1)^ ~ A ^( r ))' ( 5 - 5 ) 

CT 

n=l 

for some r/0 and for all z G C and all large t. The functions in (15.51) are polynomials in t, 
then identities (15 ,5p hold true for all t G M and in particular, at t — 1, which yields (|5.4I) . ■ 
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